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Thermalization  of  a  Fast  Ion  In  a  Plasma 
Herbert  C.  Kranzer 

PREFACE 

The  question  of  the  precise  rate  of  thermalization  of 
a  fast  ion  in  a  plasma  is  of  some  interest  In  connection 
with  several  existing  or  proposed  controlled  thennonuclear 
devices.  In  this  report,  we  follow  a  fast  ion  which  is 
injected  into  a  plasma  in  equilibrium.   Specifically,  we 
determine  the  time  history  of  the  probability  distribution 
of  this  ion  in  velocity  space.   This  is  done  by  numerical 
integration  of  the  linearized,  space-independent  Fokker- 
Planck  equation  with  both  the  ion-ion  and  ion-electron 
interaction  terms  retained.  The  mean  time  of  thermalization 
Is  calculated  for  several  widely  separated  injection  veloci- 
ties. Some  other  properties  of  this  single-ion  probability 
distribution  are  analyzed. 

The  author  is  deeply  Indebted  to  H.  Gi?ad  for  initiating 
his  Interest  in  this  problem  and  for  many  helpful  suggestions 
along  the  way.  He  is  grateful  also  for  valuable  discussions 
with  R.  Rlchtmyer  and  E.  Isaacson  concerning  numerical  pro- 
cedures, and  for  the  cooperation  of  M.  Goldstein  in  making 
available  the  necessary  computing  time. 


This  work  was  presented  at  an  AEC  -  sponsored  meeting  on 
controlled  thermonuclear  fusion  at  Gatlinburg,  Tennessee, 
April  27-28,  1959. 
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1.  The  Physical  Problem  and  Its  Differential  Equations. 

Suppose  a  single  Ion  of  velocity  ^  Is  Injected  at  time 
t  ■=  0  Into  a  homogeneous  plasma  In  thermal  equilibrium  with 
no  external  electromagnetic  fields  present.  The  probability 
distribution  f(^,  t)  of  this  Ion  In  velocity  space  satisfies 
the  space-Independent  Fokker-Planck  equation 


,2 
'r   -     --r"^s 


(^•1)       If"-^(v)-^(|^s^) 


and  the  Initial  condition 


(1.2)       t{i,0)   =  B(e  -  e^). 


(The  convention  of  sximmatlon  over  repeated  Indices  Is  employed.  ) 

The  values  of  the  dynamical  friction  coefficient  a  and 
the  dispersion  coefficient  b   to  be  taken  In  (1.1)  are  those 
corresponding  to  the  (Maxwelllan)  distributions  of  Ions  and 

ai 

3 


2 

electrons  In  the  original  plasma.   This  linearization  is 


rigorously  correct  for  a  single  Injected  Ion. 


1.  See  Rosenbluth,  MacDonald  and  Judd,  Fokker-Planck  Equation 
for  an  Inverse-Square  Force ,  Phys.  Rev.,  107,  1  (1957),  or 
Grad,  Thermonuclear  Reaction  Rates  In  an  Electrical  Discharge, 
NYO - 7977,  Inst,  of  Math.  Sciences,  N.  Y.  Univ.,  Jan.  1958. 

2.  See  Grad,  op.clt. ,  for  the  precise  form  of  these 
coefficients  In  this  case. 

3.  One  can  also  consider  f  to  describe  the  behavior  of  a 
"marked"  Ion  In  a  plasma  In  equlllbrlvun.   In  the  more  usual 
case  where  a  whole  group  of  Ions  Is  injected  or  "marked",  our 
linearization  amounts  to  ignoring  the  interaction  between  ions 
of  this  group  as  compared  to  the  forces  exerted  by  the  plasma 
on  the  group  as  a  whole.  If  one  wants  to  include  this  mutual 
interaction,  one  has  to  take  the  values  of  a^  and  bpg  corresponding 
to  a  linear  perturbation  of  the  Maxwelllan  distribution. 
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We  Introduce  the  dlraenslonless  independent  variables 
(lO)      X  =  t/(RT)^/2  , 

where  T  is  the  plasma  temperature  (assiomed  the  same  for  electrons 

and  ions),  M  is  the  ion  mass,  k  is  Boltzmann's  constant, 

R  =  k/M,  e  is  the  electronic  charge,  n  the  number  density 

of  ions  or  electrons,  and  j 

(1.5)  A=|lMl£^ 

/to  e 

is  the  ratio  of  the  Debye  length  to  the  mean  distance  of 
closest  approach  in  a  Coulomb  encounter.   (All  quantities  are 
expressed  in  cgs  electrostatic  units. )  Then  equations  (1.1) 
and  (1.2)  become  conditions 


(1.7)  0i.i,O)   =  8(5t  -  ?„)  ,    %  '  t/(M)^/2  _ 
on  the  dimensionless  probability  distribution 

(1.8)  ^{t,-,)    =  (RT)^/^  f(?,t). 


2 
Denoting  by  p  the  ratio  of  ion  mass  to  electron  mass,  the 

coefficients  in  (1.6)  are  given  (for  singly  charged  ions)  by 
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(1.9) 


with 


r    r    r  ' 


^rs   ^rs  "^  ^rs  ' 


r 


(1.10)    < 


+ 


-2^^F(x), 


L 


«r  = 


_^(p.p-i)!|p(|). 


(1.11)        p 


i  XX, 


rs 


Tl  ^  +  7;  (6 


+ 

2  ^'rs 


X  x^ 

r  s 


)  , 


where 


and 


+ 
"^1 


(1.12)    < 


^  .  %P(x), 


/?   2 


id  'dp    „/Xx 


+ 
'^2 


(1.15)    < 


=  I  tF(x)  .  ^]  , 

=  /it^^(f)^iF'(f)] 


V. 


(1.14)      P(x)  -i 


1  r  e-yV^dy 


-  e 


-xV2 


0\ir  ultimate  goal  is  to  numerically  integrate  (1.6),  (I.7), 
a  system  Involving  two  velocity  dimensions  (there  is  axial 


sywaetry  about  the  direction  of  x  )  and  one  time  dimension. 
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In  this  report  we  assvime  spherical  symmetry  In  velocity  space, 
thus  reducing  the  number  of  dimensions  by  one  and  rendering 
the  problem  more  amenable  to  nvunerical  treatment  by  a  com- 
puter of  finite  size.   If  we  are  interested  in  following  only 
the  speed  of  the  injected  test  ion,  this  should  be  an  excel- 
lent approximation. 

Hence  we  introduce  the  spherical  mean  value 

(1.15)      his  ^iV'-^)^^  =  g(x,T) 
""  |y|=x 

and  obtain  for  g  the  differential  equation 
and  the  initial  condition 


(1.17)       g(x,0)  =  x^-2  6(x  -  x^) 


The  function  Q  =  G(x)  is  defined  as 


(1.18)       G(x)  =  F(x)  +  pF(|). 
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2.   The  Numerical  Scheme 

We  proceed  to  the  ntomerlcal  solution  of  (1.16),  (1.17) 
by  finite  differences.  We  approximate  (1.16)  by  an  explicit 
difference  equation;  i.e.,  we  take  a  rectangular  mesh  with  a 
spacing  Ax  in  velocity  and  At  in  time  and  replace  ^  by  a 
forward  difference.  We  use  the  same  Ax  throughout  our  mesh, 
even  though  the  coefficients  vary  much  more  slowly  when  x  is 
large,  since  the  particular  solutions  we  seek  may  vary  rapidly 
with  velocity  for  large  x,  particularly  at  early  times. 

To  reduce  truncation  error,  the  right-hand  side  of  (l.l6) 
is  first  expanded  into  the  form 


(2.1)        |£  =  yl  .  i^  [(xG'+  G)g  +fG'  +  xG)  If  +  G  ^  ]. 


C(x)g  -f  A(x)  If  +  B(x)  ^  . 


The  X  -  derivatives  of  g  in  (2.1)  are  then  replaced  by 
centered  differences,  while  the  coefficients  A,  B,  C  are 
evaluated  analytically.  Because  A(x)  becomes  infinite  at 
X  =  0,  the  first  velocity  mesh  point  is  taken  at  x  «  -^  Ax; 
differences  centered  at  this  point  are  computed  by  assxmiing 

g  to  be  an  even  function  of  x.  The  truncation  error  of  this 

2        4 
scheme  is  of  the  order  of  (At)   or  (Ax)  . 

_p 
Since  the  maximvim  of  x  G  occurs  at  x  «  0  and  is  equal 

to  TT*  the  Courant-Friederichs-Lewy  stability  criterion 

2 
(based  on  the  ^-%  terra  alone)  is  satisfied  if 
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(2.2)        X  =  -^^  <  I/I  ^1.88. 

(Ax)"^     '^  I  "^ 

For  safety,  we  choose  X  =  1.75.  The  full  von  Neumann 
stability  criterion  is  then  satisfied  at  all  velocity  points 
except  the  single  point  x  =  -^  Ax.  This  exception  can  be 
tolerated,  since  an  oscillatory  instability  must  have  a  wave 
length  embracing  at  least  two  mesh  points.  We  thus  expect 
the  solution  of  the  difference  equation  to  be  stable  as  At 
approaches  zero. 

Incidentally,  this  does  not  mean  that  numerical  oscilla- 
tions cannot  develop.   In  fact,  the  use  of  a  centered 
difference  for  4^  will  cause  oscillations  whenever  the  effect 


of  the  "drift"  term  A  §^  in  one  time  cycle  dominates  that 

ox   2 

of  the  "diffusion"  terra  B  2_£.^  i.e.,  whenever 

Sx"^ 

(2.3)        2A^At  >  B. 


For  reasonable  values  of  At,  inequality  (2. 5j holds  when  x  is 
of  the  order  of  p  (test  ion  speed  ~  mean  plasma  electron 
speed).   In  the  one  case  that  we  consider  in  which  such 
large  values  of  x  enter,  oscillations  are  indeed  observed 
at  early  times.   What  the  von  Neumann  condition  does  insure 
is  that  any  oscillations  which  do  develop  damp  out  as  t 
increases.  This  is  exactly  what  does  happen. 

In  most  of  the  cases  considered  below,  the  initial  values 
(1.17)  of  g  are  represented  numerically  by  choosing  x^  to 
coincide  with  a  mesh  point  and  taking  g(x,0)  to  be  x  ~  (Ax)~ 
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at  this  mesh  point  and  zero  at  the  others.   In  the  one 
case  where  x  >  p,  this  prescription  would  cause  the 
oscillations  mentioned  in  the  preceedlng  paragraph  to  reach 
unacceptably  large  amplitudes  before  damping  out.   In  this 
case,  we  begin  the  calculation  at  a  time  t  >  0  which  is 
small  relative  to  the  time  scale  of  thermalization.  The 
starting  values  of  g  are  chosen  to  be  the  values  g  would 
have  at  time  t  if  the  coefficients  A,  B,  C  in  (2.1)  had 
everywhere  the  values  they  have  at  x  : 

exp  [Ct  -(x-x  +  At)^/4Bt] 

(2.M        g(x,^)  =  .-^ . 

2x^  (ttBt)  ' 

The  choice  of  numerical  upper  bounds  for  x  and  t  is 
facilitated  by  the  fact  that  the  exact  solution  of  (1.16), 
(1.17)  approaches  as  t  — >  «  the  Maxwellian  distribution 


(2.5)        gw(x)  =/|e-^^/2 


Thus  for  an  x  not  much  larger  than  x  the  solution  remains 
altogether  negligible  for  all  time.  We  take  such  an  x  for 
an  outer  boundary  and  impose  there  the  artificial  boundary 
condition  g  =  0.  The  upper  bound  on  t  is  determined  by  the 
computation  itself:  we  stop  computing  whenever  g(x,T)  is 
as  close  to  giy,(x)  as  the  numerical  approximation  allows. 
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J>.      Injection  Velocities  and  other  Parameters. 

Nvimerlcal  computations  have  been  performed  for  a  plasma 
consisting  of  deuterium,  for  which  p  =  60.5948.  Four  widely- 
spaced  Injection  velocities  x  were  chosen:  x  =  0, 

o  o 

Xq  =  1.55,  Xq  =  9-7,  and  x^  =  240.  The  first  of  these  was 
chosen  to  provide  a  reference  relaxation  time.  The  second 
corresponds  to  the  center  of  the  Initial  distribution  used 
by  Rosenbluth's  group  In  a  Fokker-Planck  equation  calcula- 
tion of  the  thermallzatlon  of  a  group  of  test  Ions  In  the 
absence  of  an  external  plasma.  The  choice  of  the  last  two 
Injection  velocities  can  be  most  easily  understood  by  reference 
to  the  "friction  curve"  --  the  plot  of  the  total  friction 
coefficient  a  (cf.  (1.6))  as  a  function  of  x  (see  Fig.  1). 
This  curve  begins  at  zero  (for  x  =  0),  rises  to  a  maximum 
value  Ow  at  about  x  =  1.3,  decreases  to  a  minimum  at  x  =  9.7, 
rises  to  a  second  maximum  value  -p  <%  at  x  =  1.5p  '^  80,  and 
finally  falls  off  toward  zero  In  proportion  to  x  .   The 
velocity  X  =  9«7  lies  at  the  central  minimum,  while  x  =  240 
Is  approximately  that  velocity  for  which  a  reaches  on  Its 
final  downward  curve  the  value  that  it  has  at  x  «  9.7. 

The  velocity  mesh  width  Ax  must  be  taken  small  enough 
to  yield  a  significant  nxomber  of  points  in  the  region 
0  <  X  <  1.3  in  which  a  rises  to  its  first  maximum.  This  is 
also  the  region  in  which  coefficients  p__  (particularly  p^„ ) , 


4.  MacDonald,  Rosenbluth,  and  Chuck,  Relaxation  of  a  System 
of  Particles  with  Coulomb  Interactions,  Phys.  Rev.,  I07,  35O 
(1$57). 


-li- 
as well  as  the  limiting  solution  g]^(x),  are  substantially 
different  from  zero.  Therefore,  Ax  was  never  chosen  larger 
than  0.1. 

The  case  x  =  240  requires  smoothing  of  the  initial 
data,  as  described  in  Section  2.  The  time  ^   of  smoothing 
(see  (2.4))  was  taken  as  37-7. 

The  nximerical  parameters  used  in  the  four  cases  are 
summarized  in  Table  I.  The  computations  were  performed  on 
the  IBM  704  at  New  York  University. 


Injection  velocity  x 

0 

1.55 

9.7 

240 

Initial  smoothing  time  t 

0 

0 

0 

37.7 

Velocity  mesh  ZIvx 

0.1 

0.05 

0.06 

0.1 

Time  mesh  At 

0.0175 

0.004375 

0.0063 

0.0175 

Upper  velocity  bound  x 

7.0 

7.0 

12.0 

250 

Table  I  —  Details  of  the  numerical  computation 
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4.  The  Numerical  Results. 

In  the  three  cases  where  x  is  positive,  the  solutions 
behave  quite  similarly  (see  Fig.  2).  At  first,  the  Initial 
delta -function  diffuses  into  an  asymmetric  Gaussian  distribu- 
tion of  varying  width,  steeper  on  the  high-velocity  aide 
(Fig.  2A).  The  peak  of  this  distribution  migrates  toward 
lower  velocities  (Fig.  2B).  This  regular  evolution  continues 
until  the  inner  tall  of  the  Gaussian  reaches  x  =  0  (Fig.  2C). 
At  that  time,  a  second  peak  forms  at  x  =  0  (Fig.  2D).  This 
peak  rapidly  Increases  in  height  and  width  (Fig.  2E),  completely 
swallowing  up  the  original  Gaussian,  until  (Fig  2F)  it  some- 
what surpasses  the  Maxwelllan  distribution  (2.5).   Finally 
(Pig.  2G)  it  slowly  falls  back  and  broadens,  approaching  (2.5) 
asymptotically. 

When  X  =  0,  the  initial  distribution  is  Just  an  extreme 
case  of  the  pattern  of  Flg\ire  2F.  Hence  the  asymptotic  decay 
toward  the  Maxwelllan  distribution  begins  immediately. 

We  proceed  to  give  quantitative  details  of  this  general 
picture  in  the  various  cases.   A  physical  idea  of  the  size 
of  the  units  involved  may  be  obtained  by  noting  that  the  root 
mean  square  plasma  ion  speed  corresponds  to  x  =  V^-^l.T^* 

while  the  r.m.s.  electron  speed  corresponds  to  x  =p /3  <»>  105. 

5 
Furthermore,  the  Spitzer  ion-ion  collision  tlme-^ 


(4.1)       t^  -  M^/^(5kT)^^V  8  X  0.714  Tme^lnA 
corresponds  to  t  »  3.62. 


5.   Spitzer,  Physics  of  Fully  Ionized  Gases,  Intersclence, 
New  York,  1955- 
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The  early  stage  of  thermalizatlon  (through  Fig.  2C) 
may  best  be  understood  by  reference  to  the  solution  (2.4) 
of  the  equation  obtained  from  eq.  (2.1)  by  making  Its 
coefficients  constant.  This  solution  represents  an  exact 
(symmetric) Gaussian  distribution  whose  width  Is  proportional 
to  /Bt  and  whose  peak  moves  toward  lower  values  of  x  with 
the  constant  speed  A.   Since  the  actual  coefficients  in 
(2.1)  vary  slowly  with  x,  at  least  for  x  >  1.3,  it  is  not 
surprising  that  the  correct  solution  turns  out  to  be 
approximately  Gaussian  at  early  times.   Indeed,  the  motion 
of  the  peak  toward  lower  x  takes  place  with  an  instantaneous 
speed  which  is  approximately  equal  to  the  value  of  A  at  the 
Instantaneous  position  of  this  peak.   (See  Fig.  J>,   where 
the  maxima  and  minima  of  the  friction  curve  are  clearly 
mirrored  in  the  slopes  of  the  peak-x-versus-x  curves.) 

Even  the  slight  asymmetry  of  the  solution  can  be  explained 
on  this  basis.   One  need  only  observe  that  B(x)  is  a  monotone 
decreasing  function  of  x,  so  that  the  distribution  should 
have  a  greater  "width"  to  the  left  of  its  peak  than  to  the 
right . 

In  Figure  4,  we  plot  the  width  (distance  between  half- 
peak  points)  of  the  solution  as  a  function  of  the  position 
of  the  peak  in  two  of  our  cases.   Note  that  the  width  does 
not  depend  simply  on  the  monotone  dispersion  p,  but  rather 
seems  to  vary  with  the  friction  coefficient  (or  perhaps  with 
dp/dx).   Note  furthermore  that  in  both  cases  shown,  especially 
in  the  case  x  =  240,  the  peak  width  in  the  early  stages 
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surpasses  the  ultimate  width  2  /2/n2  »  2.4  of  the  limiting 
Maxwelllan  distribution  g^Cx). 

The  graphs  In  Figures  3  and  4  are  changed  to  dotted 
lines  at  the  time  t  =  x^  when  the  central  peak  of  Pig.  2D 
begins  to  form.  They  are  cut  off  completely  at  t  =  t^,  the 
time  at  which  the  original  peak  no  longer  exists  (between 
Figs.  2E  and  2F).   We  call  period  t,  <  t  <  Tp,  during  which 
the  central  portion  of  the  Maxwelllan  distribution  (2.5) 
is  filled  out,  the  middle  stage  of  thermalization.   In 
Figure  5  we  show  the  heights  of  the  central  peaks  as  a  func- 
tion of  time  during  this  middle  stage  and  continuing  on  to 
(and  in  some  cases  beyond)  the  time  t  =  tl  defined  below. 

The  late  stage  of  thermalization,  beginning  at  t  =  Tp, 
is  characterized  by  the  asymptotic  decay  of  the  solution 
toward  the  Maxwelllan  distribution  gw(x).  From  physical 
reasoning  (cf.  Spitzer-^),  it  is  plausible  to  expect  that  this 
decay  becomes  exponential  at  late  times.   Accordingly,  we 
compute  numerically  the  average  deviation 

00 

(4.2)        6(t)  =  I  x2|g(x,T)  -  gp,(x)|  dx 

o 

of  the  solution  from  Maxwelllan  and  plot  this  deviation 
against  time  on  a  semllogarithmlc  graph  (Figure  6).  We 
observe  that,  in  all  four  cases,  this  graph  approximates  a 
straight  line  from  some  time  t  =  t,  >  Tp  on  to  quite  late 
times.  Thus  we  can  conclude  that  there  exists  in  each  case 
a  "relaxation  time"  t  such  that 
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(4.3)       C(t)  =  const,  e"*^^"^  for  t  >  x^. 


(The  aberrations  in  the  graphs  of  Figures  6a  and  6b  at 
very  late  times  are  due  to  purely  numerical  causes:  6(t) 
becomes  less  than  the  truncation  error.  The  leveling  off 
In  Figures  6C  and  6d  represents  the  filling  out  of  the 
Maxwell  tall.  ) 

We  define  the  mean  time  Tj^  of^  thermalizatlon  of  an  ion 
to  be  the  time  at  which  the  average  deviation  of  g(x,T) 
from  a  Maxwellian  distribution  is  20  percent  (i.e.,  d  =  0.2). 
The  times  -x-^,   Tp,  t,,  t^,  t  are  summarized  in  Table  II.   In 
the  case  x  =  240,  the  smoothing  time  t  is  Included  in  all 
elapsed  times. 

Also  shown  in  Table  II  is  the  maximum  deviation,  due 
to  truncation  and  roundoff  errors,  of  the  total  nxamber  of 
test  particles 


(H.H) 


x  g(x,T)  dx 


J 

o 


from  its  theoretical  value  1  during  the  time  0  <  t  <  -X).. 
The  small  size  of  these  deviations  is  an  indication  of  the 
accuracy  of  the  computation. 
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Injection  velocity  x 

0 

1.55 

9.7 

240 

Time  T,,  at  which  central 
peak  begins  to  form 

0.9 

50 

1840 

Position  X  of  outer  peak 
at  time  x. 

1.07 

7.08 

11.75 

Time  Xp,  at  which  outer 
peak  disappears 

1.3 

100 

1970 

Position  X  of  outer  peak 
at  time  t^ 

0.82 

3.90 

4.2 

Height  g  of  outer  peak 
at  time  t^ 

0.55 

0.013 

0.0085 

Time  T-,,  at  which  exponen- 
tial decay  begins 

0 

3.9 

150 

1990 

Relaxation  time  x  of 
exponential  decay 

3.4 

6.7 

16 

23 

Mean  thermalization 
time  T2^ 

6.0 

4.4 

157 

2043 

Ratio  of  x^  to  Spitzer 
collision  time 

-2 

-1.2 

^45 

-575 

Maximum  deviation  of  particle 
number  from  unity  (up  to 
time  Xi^) 

.001 

10"5 

.001 

.004 

Table  II  -  Principal  numerical  results 
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5.   Conclusions. 

The  thermallzatlon  times  given  In  Table  II  were  the 
main  goal  of  this  Investigation.   Several  other  points  also 
deserve  mention,  however. 

First,  some  trivial  observations.   It  seems  clear  from 
the  slope  of  the  curve  in  Fig.  5A  that  an  ion  whose  speed 
X  »  9-7  Is  slowed  down  principally  by  encounters  with 
electrons,  while  for  x  «  9.7  it  is  the  other  ions  that  do 
the  Job.   For  x  fni'^.'J,   both  are  important.   Furthermore, 
Fig.  4  shows  that  an  ion  is  extremely  unlikely  to  gain  energy 
unless  it  is  already  below  the  mean  energy  of  the  plasma 
particles. 

Secondly,  no  surprise  should  be  occasioned  by  the  dif- 
ference in  the  relaxation  times  t  in  the  various  cases.  The 
relaxation  time  for  x  =  0  agrees  fairly  well  with  the 
Spitzer  collision  time.  The  other  relaxation  times  are 
longer  because  the  corresponding  distributions  are  broader 
when  they  start  to  relax  (see  line  k   of  Table  II).  Thus  the 
Spitzer  collision  time  cannot  be  taken  as  an  absolute  measure 
of  the  rate  of  approach  to  equilibrium. 

Of  course,  the  existence  of  an  exponential  decay  toward 
a  Maxwellian  distribution  does  not  contradict  Rosenbluth's 
observation  that  the  Maxwell  tail  takes  infinitely  long  to 
fill  up.  The  reason,  as  (4.2)  shows,  is  that  the  tall  makes 
a  contribution  to  the  deviation  <3(t)  which  decreases  faster 


than  exponentially.   In  fact,  the  point  x  =  X  of  the  tall  is 

I2     ?  -y  /p 
weighted  with  a  factor  y—  X  e  ^  '  "^     in  the  computation  of  6. 
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Finally,  It  Is  instructive  to  compare  our  niimerlcal 
solution  in  the  case  x  =  1.55  with  the  solution  of  Rosenbluth 
mentioned  in  Section  3  above.  This  comparison,  made  in  Table 
III,  shows  a  qualitative  similarity  in  the  mode  of  thermallza- 
tlon  in  these  two  extremely  different  situations.  The  speed 
of  therraallzation,  however,  is  almost  twice  as  great  in  the 
present  calculation.   In  other  words,  if  a  group  of  ions, 
originally  all  at  the  same  speed  i,   come  to  Maxwellian 
equilibrium  under  the  forces  between  them  in  time  T,  then  an 
ion  injected  with  velocity  i   into  the  resulting  plasma  will 
thermallze  (i.e.,  attain  a  Maxwellian  velocity  distribution) 
in  a  time  slightly  greater  than  t/2. 


Source 

This  report, 
^o  =  1-55 

MacDonald, 
Rosenbluth, 
and  Chuck^ 

External  plasma  present 

yes 

no 

-^1 

0.9 

1.5 

X  of  peak  at  t. 

1.07 

1.03 

-^2 

1.5 

1.75 

T 

6.7 

11 
(estimated) 

Table  III  -  Comparison  of  time  scales  of 
thermallzatlon  in  the  presence  and  absence 
of  an  external  plasma. 
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